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INVESTIGATION OF THE DYNAMIC RESPONSE 
OF AIRPLANE WINGS TO GUSTS 
By Harold B. Pierce 


SUMMARY 


A method of predicting the dynamic response of airplane wings 
to gusts "by considering only the fundamental mode of tending is 
presented, together with the results of model tests made to evaluate 
the method. In addition, the results of a series of calculations 
obtained "by using the method are given to illustrate changes in 
the dynamic response of ai'rplane wings "brought about by changes in 
gust and airplane parameters . An appendix giving the details of the 
method and the procedure for rh® determination of constants is 
included. 

Although the te3t results are not suitable for predicting 
dynamic-stress ratios, they serve to indicate that the method is 
of sufficient accuracy to predict the ratio of maximum dynamic 
wing deflection to maximum fuselage acceleration increment for 
conventional airplanes. The test results also illustrate the 
need for including aerodynamic damping in calculations of dynamic 
response of airplane wings. 

The calculations made to show the effect of change of certain 
gust parameters indicate that: 

(1) The dynamic -stress ratio for airplane wings increases as 
the gradient distance of the gust decreases. 

(2) For the assumed design gust of 10-chord gradient distance, 
the overstress in a single gust may be as much as 12 percent. 

( 3 ) Although the results for one type of repeated gust are 
not conclusive, a repeated gu3t does not 3eem to be more critical 
than a single gust. 
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The calculations made to show the effect of change of certain 
airplane parameters indicate that for the assumed design gust 
of 10-chord gradient distance: 

(1) A change in forward velocity of the airplane does not 
appreciably change the dynamic -stress ratio. 

(2) A reduction of wing frequency by a change either in weight 
or in stiffness of the wing results in an increase in the dynamic - 
stress ratio. 

(3) The wing-tip acceleration increment is generally much 
greater than the fuselage acceleration increment and the ratio of 
the two tends to increase as the 3 peed. increases. 


INTRODUCTION 


In the present-day design of airplane wings for strength, the 
predominating loads considered are caused by maneuvering and by 
gusts. For transport-type airplanes, for which the design maneu- 
vering loads are relatively low, the design gust loads are fre- 
quently critical, especially as the operating speed increases. 

These gust loads, for- simplicity of calculation, are assumed to 
be static loads applied to a rigid airplane. For the smaller 
airplanes, the assumption of a rigid structure appears to be 
reasonable. With increase in size, however, airplane wing3 tend 
to deflect more and, thus, the assumption of rigidity is subject 
to question. Consequently, the trend toward larger and faster 
airplanes necessitates the determination of methods for calcu- 
lating and evaluating the effects of dynamic response. 

In the pa 3 t a number of investigations of the dynamic response 
of airplane wings on encountering gusts have been conducted. 

Notable among these investigations are the ones reported in 
references 1 to A. In reference 1, under the assumption of no 
pitching motion, Kussner sets forth the exact equations for the 
loads and moments when the airplane encounters a gust and applies 
these equations to the flexural system which is composed of the 
wing and fuselage of en airplane . Since a mathematical error was 
found in his simplification of these integro-differential equations, 
the calculated results given by Kussner are therefore considered 
invalid. References 2 and 3 present two other simplifications of 
the problem, but the dynamic -stress results are difficult to 
dissociate from the stability reactions in which they are contained. 
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Reference b presents still another simplification that is "based on 
an infinite-mass fuselage. These reference papers represent 
isolated investigations of variov.s parameters affecting the 
dynamic response of airplanes to gusts and use a number of methods 
that cannot readily be evaluated by comparison « In addition, the 
validity of the various investigations has never been experimentally 
shown . 

Upon consideration of the need for investigating the flexural 
reactions of large airplanes in gusts and because of the omissions 
in the existing literature, an investigation was undertaken with 
the following purposes in mind: 

(1) To evolve a relatively simple, yet sufficiently accurate 
method for determining the dynamic stress in the wing3 of a given 
airplane due to their dynamic response tc gusts. 

(2) To make an experimental check of the method. 

(3) To make a study of the effect of the changes of gust and 
airplane parameters on the stress due to the dynamic response of 
airplane wings . 

The results of this investigation are presented herein, together 
with the analytical method, which is restricted to the fundamental 
mode of bending and includes aerodynamic damping. 


SYMBOLS 




equivalent mass of wing, slugs 


e 


Mf 


equivalent mas3 of wing-fuselage, slugs 


e 


M 


mas3 of airplane, 3lugs 



absolute displacement of equivalent wing, feet 


5 f 


absolute displacement of equivalent wing -fuselage, feet 


wing-tip deflection increment, feet, (By - Bf) 
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^d-max 

8 st 




X 


K 


f 

t 


F(t) 

F„ e (t) 

F f (t) 
e 


Ally 


An T 


max 


W 

W 

w e 

U 

Umax av 

AO 


maximum value of (5 W - 6f) , feet 

deflection of equivalent ving under conditions corre- 
sponding to normal static design procedure, feet 

equivalent ving damping coefficient, pound-seconds per 
foot 

equivalent ving-fuselage damping coefficient, pound-seconds 
per foot 

damping coefficient of entire airplane, pound-seconds 
per foot, ^ X w ^ + Xf^ j 

equivalent spring constant "based on ving frequency and 
equivalent ving mass, pounds per foot 

ving frequency, cycles per second 

time, seconds 

forcing function on entire airplane, pounds 

forcing function on equivalent ving, pounds 

forcing function on equivalent ving-fuselage, pounds 

acceleration increment normal to chord of ving, g uni ts 

acceleration increment on rigid airplane, g units 

maximum value of acceleration increment on rigid 
airplane, g units 

proportion of total air load assumed to deflect equivalent 
ving 

veight of airplane and of equivalent "biplane, pounds 

weight of equivalent ving, pounds 

gust velocity, feet per second 

average maximum gust velocity, feet per second 

pitch -angle increment, degrees 
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V forward velocity, feet per second 

D differential operator 

Hp gradient distance of the first gust (fig. 2), chords 

Hp gradient distance of the second gust (fig. 2) , chords 

2-^ distance "between the end of the gradient distance of first 

gust and the "beginning of the gradient distance of 
second gust (fig. 2), chords 


An& W 
te- bt 


at 


"bt 


1 


= VbeAria 


An a arbitrary load-factor increment that the airplane would 

experience if if had no vertical motion when traversing 
the gust, g units 

"b = ^ at bt = 1 or at the maximum value of the function te'^ 


ANALYTICAL METHOD 


Examination of the problem of the dynamic response of airplanes 
to .gusts shows that the spanwise distribution of the imposed loads 
varies over wide limits (reference 5) and that the structure itself 
may react in a number of combinations of' modes. "Rather than attempt 
to 3olve the general case, consideration was given to a method for 
determining the origin of the principal stresses with the intention 
of reducing the problem to one of reasonable dimensions. A3 a 
result, this paper is concerned with the analysis of the fundamental 
mode of bending of the wing under symmetrical loads which are 
assumed to be known. 


Development of Method 

Eu.ie_i.lyj tue present method follows the steps given in 
references 2 and 3 in that the airplane is replaced by a simple 
equivalent aerodynamic and elastic syst- i. The "basic forms of 
forcing iunct ion to he used for the response of an airplane to 
known single gusts, however, have been selected from the results 
of gust-tunnel tests such as those described in reference 6, 
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Biplane as equival e nt airplane .- The problem of determining 
the stress caused by dynamic response is simplified by the previous 
assumption that the fundamental mode of vibration of the ving is 
the most significant. As a result, the deflection of the ving tip 
with respect to the fuselage may be taken as a direct measure of 
the stress in the wing. The equivalent system, then, is represented 
as a biplane (fig. 1) that has the motions of the rigid upper wing 
and the rigid lower wing-fuselage combination adjusted to have the 
same motions as the ving tip and fuselage of the airplane under 
investigation. The equivalent biplane system must include the 
proper distribution of aerodynamic as veil as inertia and elastic 
forces or force coefficients. The problem is then resolved into 
one of obtaining the proper constants to be used as coefficients 
in simultaneous linear differential equations which represent the 
equivalent biplane. These equations follow: 


d6 w 

Mv e + H d T + K(5v ’ 6f) “ Fw e (t) (1) 

CL o 


d c 5,» d5.p 

+ e IT * K(5 v ' 5 f> = (2) 


These equations are solved for the stress in the wing as represented 
by the wing deflection (5 W - 5 f ) and for the load parameters 
represented by the fuselage and wing-tip acceleration increments. 

The details of the solution and the procedure for determination of 
the coefficients ar© given in the appendix. 

Since the normal gust-load design procedure assumes that the 
lead is applied statically and that the normal acceleration is 
constant along the span of the airplane, the dynamic wing deflection 
must be compared with the static wing deflection under the action of 
the acceleration of the "rigid" airplane. In order that the comparison 
be valid, the rigid-airplane acceleration iCn r is determined as 
the response of the equivalent biplane to the sum of the forcing 
functions F v e (t) and Ff e (t) when the springs are replaced by 
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rigid rods. The static deflection 5^ is then calculated from 
the simple formula 


which, in accordance with normal design procedure, considers no 
aerodynamic damping of the static wing deflection. The stress 
computed hy the static-load method is therefore multiplied "by the 


Of further interest is the fact that, since equations 1 and 2 
are linear, ratios of the maximum flexible-wing acceleration 
increments to the corresponding maximum rigid-airplane acceleration 
increments are independent of the magnitude of the load and may "be 
applied directly as multiplying factors to those acceleration 
increments computed in the determination of the gust load. 

Required conditions .- In arriving at the characteristic values 
of the coefficients for the equivalent biplane, the following 
conditions should be satisfied: 

(1) The total mass of and the total load on the equivalent 
biplane should be identical with those of the original airplane. 

(2) The kinetic energy of vibration of the upper wing should 
closely approximate that of the original wing beam for an amplitude 
of vibration of the upper wing equal to that of the tip of the 
original wing. 

(3) The natural frequency should be the same as that for the 
fundamental mode of the original wing. 

(4) The upper wing should deflect under the equivalent static 
load the same amount that the wing tip of the airpiane does under 
the corresponding aerodynamic static-load distribution. 

(5) The damping coefficient of the upper wing should represent, 
at least to peak load, the damping of the motion of the wing. 


Eq uivalent masses and spring constant .- The equivalent masses 
and spring constant are determined from conditions (l) to (3) . 

The equivalent mass of the wing is generally determined 

G 



(3) 


ratio 5 dma3C /5 e f 


Discussion of Terms 
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firat and the equivalent spring constant K is then determined 
from the known or estimated wing frequency. The equivalent 
ving-fuselage mass Mf is then simply the total mass of the 

airplane M minus M Wp . For special cases, such as strut- 

supported wings and arrangements other than the conventional 
cantilever monoplane, considerable care is required. 

Damping coefficient ..- The damping coefficient is actually 
the lift per unit vertical velocity for each wing section along 
the span of the airplane and, as such, includes the effect of the 
vertical motion of the airplane as a whole, as well a3 that of 
the vibratory velocity caused by the bending action of the wing. 
Although the lift or damping force arising from the two motions 
would be subject to the effects of unsteady flow, an analysis 
and some tests, described in the appendix, indicate that an 
acceptable approximation to the unsteady-lift damping load is to 
represent this load as 75 percent of the steady-lift damping load. 
The proportion of the total damping coefficient assigned to the 
upper wing of the equivalent biplane should change as the airplane 
traverses the gust because of the change in the relative signifi- 
cance of the wing-deflection velocity and the vertical velocity of 
the airplane as a whole. In a given calculation, however, the 
proportion used is assumed to be determined by the relation 
existing between the velocities at the time of maximum wing 
deflection. Since this relation is not known prior to the calcu- 
lations, an approximate criterion, described in the appendix, 
was determined for the division of the total damping coefficient X 
into X-y^ and Xp « 

Forcing functions .- The forcing function is considered to be 
the time history of the air loads applied to the rigid airplane 
minus the damping load due to vertical motion. The omission of 
the damping load due to vertical motion from the total air loads 
is necessary because it is included in the damping terms of 
equations (l) and (2), which provide for approximation of both the 
vertical motion and vibratory damping. With this omission in mind, 
the forcing function may be determined by any of several available 
methods. One method used by Kussner (reference 1) is to set up the 
basic equations in terms of gust velocity and. its spanwise distri- 
bution and to include these equations directly as the forcing 
functions in equations (l) and (2) . A second method would he to 
solve a set of equations such as given in reference 7, which 
describe the reaction of a rigid airplane to a guat and, from the 
resulting time history of acceleration, to determine a curve 
representing the forcing function F(t) . In order to permit this 
second form of the forcing function to be used readily, a solution 
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is made of equations (l) end (2) with a unit-function type of 
forcing function (reference 8) for a desired reaction such as 
wing deflection. The results are then combined, through the use 
of Duhamel's integral which may be solved graphically through 
use of the method given in reference 9* 

Since for these calculations the shape of the curve repre- 
senting the forcing function is the principal characteristic and 
since the unsteady-lift effects and the effects caused by stability 
reactions are of doubtful accuracy insofar as their predictions 
are concerned, it seemed desirable to make 3ome approximation of 
the forcing function to obviate its calculation for each airplane. 
This approximation was accomplished through the analysis. of 
accelerometer records from the teats of rigid-airplane models 
traversing known gusts of the type shown in figure 2(a) . The 
procedure followed was to compute the time history of the vertical 
velocity during the traverse of the gust, to determine the damping 
force due to this vertical velocity, and to add the acceleration 
increment due to this damping force to the net acceleration incre- 
ment from passage through the gust, For representative gust 
shapes, from a sharp-edge gust to one rising to maximum velocity 
in 20 chords, the resulting time histories of the forcing functions 
were determined and it was found that, for practical purposes, all 

curves could be represented by a function of the type Ate”^ 1 " . 

When unconventional designs are considered, however, it appears 
that the second method mentioned or tests of the model in the 
gust tunnel should be used to determine the forcing-function type , 

The specific 3hape of the curve representing the forcing 
function Ate ^ for a given calculation is determined from aa 
approximate relationship between the chordwise velocity distri- 
bution of the chosen gust and the reaction of the given airplane 
to the forcing function. The tests on rigid-airplane models also 
showed that the maximum value of accelera-tion increment occurs 
approximately at the time that the model reaches the end of the 
gust gradient distance H 1 . As a result, the choice of the 

value b in the forcing function Ate"^ to make this forcing 
function represent a given gust gradient distance is merely one 
of making the time history of the rigid-airplane acceleration 
increment An r reach its maximum value at the same time that the 
airplane reaches the end of the gust gradient distance H^. A 

general procedure for the selection of values of b to represent 
given gusts is outlined in the appendix. 
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The forcing function F w ^ ( t ) is ascertained "by determining 

the equivalent static load on the upper wing of the equivalent 
"biplane which would yield the 3 ame deflection of the upper wing 
as the total distributed static air load yields for the original 
wing tip. The forcing function F f (t) is then merely the dif- 
ference between F(t) and F w (t) . 


EXPERIMENTAL INVESTIGATION 


In order to obtain informed ion on the reliability of dynamic 
response calculations of the type considered herein and the dif- 
ficulties to be expected in any actual application, tests were 
made in the gust tunnel of a simple model equipped with semirigid 
wings. These tests were conducted for three conditions of wing 
stiffness and for three gust-gradient distances to determine the 
fuselage acceleration and the corresponding wing-tip deflection 
increments . 


Apparatus 

The model used for the tests is shown in figure 3. Pertinent 
characteristics of the model are given in table I and in figure h. 
The ving panels were rigid and connected to the cabane by flexure 
plates to permit freedom of motion in "bending" while offering a 
maximum stiffness in torsion. Struts connected by universal joints 
to the wing (fig. ?) were supported on the inboard end by canti- 
lever springs which could easily be changed to modify the na-tura.1 
wing frequency. The rigid-wing condition was obtained by attaching 
the struts directly to the fuselage. 

The model carried a miniature accelerometer and lights at the 
nose and tail and other lights were located at the wing-strut 
connections to impress records of the wing-tip deflection on the 
accelerometer film through a lens mirror system as indicated in 
figure 6. A sample accelerometer record is shown in figure 7 where 
the time histories of wing-tip deflection are labeled A and the 
acceleration time history is labeled B. The distance C between 
the time histories results from the fact that the two record light 
beams strike the film 90° apart on the drum. 

The gust tunnel and associated equipment have been described 
in detail in reference 6. The gust types pertinent to the present 
paper are shown, together with test results subsequently described, 
in figures 8 to 10. 
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Tests 

Tests consisted of a minimum of 10 flights for each of the 
wing-frequency values, 13 . 5 > 26.1, and infinite cycles per second, 
end related spring constants and for each of three gust types. The 
tests were made for one forward velocity, gust velocity, and air- 
plane weight so that the only variables were gust type and wing 
frequency. For each flight of the model through the gu3t, measure- 
ments were made of the forward velocity, gust velocity, normal 
fuselage acceleration, wing-tip deflection, and the pitching 
motion. 


Precision 


The measured quantities are estimated to be accurate within 
the folloxving limits for any flight: 


Normal fuselage acceleration increment, g units t0.05 

Wing-tip deflection, inches . +0.01 

Forward velocity, feet per second ..... +1.0 

Gust velocity, feet per second +0.1 

Pitch-angle increment, degrees +0.2 


In addition to errors in the recorded quantities resulting from 
instrument characteristics and limitations on reading accuracy, 
another error which is difficult to evaluate is that of the uni- 
formity of the gust shape both laterally and longitudinally. The 
longitudinal variation (in the direction of flight) of the gust 
shapes shown in figures 6 to 10 would tend to modify the shape of 
curve representing the forcing function from that assumed for a 
linear gradient. With such variations existing in the direction of 
flight similar variations micht exist across the span of the model. 
Inspection of more recent surveys taken 1 foot on either side of the 
center line indicate that, at points in the gust, such variations 
could be of the order of 12 percent of the center value. 


Results 


The launching is intended to set the model in a steady straight 
glide at constant speed without pitch, roll, or yaw. In actual 
tests the fulfillment of these ideal conditions is almost impossible 
anu therefore, all flights in which pitch, roll, or yaw were excessive 
prior to entry into the gust must be disregarded. The records from 
each satisfactory flight were evaluated to obtain time histories 
of acceleration increment, wing-tip deflection increment, and pitch 
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increment during the entire traverse of the gust and to obtain the 
forward speed and gust velocity for the flight. Sample time his- 
tories of the acceleration increment and the wing-tip deflection 
increment are shown in figures 8 end 9 for the two wing frequencies 
used, 13*5 and 26.1 cycles per second, respective ly, and the time 
histories of acceleration increment are shown in figure 10 for 
the rjgid-wing condition. The results are plotted against the 
distance penetrated into the gust in chords and have been corrected 
to a forward speed of 6l feet per second and a gust velocity 
of 6 feet per second on the assumption tha+ the acceleration 
increment and wing-tip deflection increment are directly propor- 
tional to forward speed and gust ■velocity for small variations in 
these quantities. Results for two flights under similar con- 
ditions are 3hown as a sample of the data obtained for each 
condition of wing stiffness and gust shape. 

In order to provide a measure of the pitching response of the 
airplane model, data are pr< sented giving the pitch-increment 

ratio at the instant of maximum acceleration increment. 

U / v 

Analysis of past gust-tunnel tests has shown that this ratio is 
approximately equal to the fraction by which the measured 
acceleration increment differs from that predicted by the analysis 
of reference 1C, which assumes that the airplane does not pitch 

while traversing a gust. Average values of for e&oh test 

condition are shown in figure 11 as a function of the distance to 
peak acceleration increment to indicate the dependence of the ratio 
on gust gradient distance. Also shown js the scatter band of data 
for four moaels (six conditions of weight end speed) from the tests 

of which the forcing function Ate""* 3 * was derived. 

Figure 12 shows the maximum values of acceleration and 
wing-tip deflection increments for each wing frequency, for all 
records read, as the wing -tip deflection per unit acceleration 
plotted against the gradient distance. The right and left wing 
results have been given separately to show the degree of symmetry 
of the wing-tip deflections. 


CALCULATIONS 


In the preceding sections of this report, a method, was outlined 
for calculating the dynamic response of airplane wings under the 
action of arbitrary gusts ar ’ some simple tests of a model in the 
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gust tunnel were described. The .remaining problems are to determine 
whether, considering the experimental data, the method yields 
results of reasonable accuracy and, finally, to investigate in a 
general way some of the significant factors which determine the 
degree of dynamic response of airplane wings in turbulent air. 

The first problem - that of checking the accuracy of the 
method against experimental data - may be solved simply by checking 
the data obtained for the airplane model against calculations 
based on the conditions for which the tests were made. Data for 
the solution of the second problem, however, requires calculations 
for representative airplanes for which dynamic response might be 
of some concern. The airplanes selected for these calculations 
were those for which some information as to structure and mass 
distribution was available. 

Briefly, the calculations were made to indicate the effects 
of airplane weight, weight distribution, gust size, airplane speed, 
and wing stiffness and to investigate the possibility of further 
simplification in the method of calculation. 


Calculations for Comparison with Experimental Results 

In the dynamic -stress calculations for the test model the 
principal differences from the method 3et forth in this paper for 
the normal cantilever-wing airplane arise from the fact that the 
wings of the model deflect about the hinge as rigid bodies. As a 
result, parameters affecting angular-frequency reactions pre- 
dominate. In other words, the actual spring constant of the model 
is in units of torque per degree of deflection and the reactions 
of the wing also depend on the moment of inertia of the wing panels 
about the hinge point. The moment of inertia of the wing panels, 
therefore, was determined experimentally and the equ.iva.lent wing 
mass for the calculations was that mass which, if placed at the 
model wing tip, would have the same moment of inertia about the 
hinge points as the wing panels . Since the wing frequencies for 
the two model conditions were known, the equivalent linear spring 
constants were then determined. 

Further differences from the normal cantilever- wing case arise 
when consideration is given to the division of the air and damping 
loads applied to the model into the corresponding loads for the 
equivalent biplane. The division of the damping coefficient in 
the sharp-edge gust and in the gust with 10-chord gradient distance 
can be shown to be approximately 0.50 and 0.50 for the equivalent 
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ving and wing-fuselage, respectively, whereas for gradient distances 
greater than 10 chords the division approximates that for the air 
load, namely, O.37 and O.63 for the wing and wing -fuselage. 


The equivalent constants, together with the values such as 
test forward velocity and the gradient distances of the test 
gusts, are included in table I. The actual gust profiles 

' U/UmaT plotted against horizontal distance from edge of gust 
\ ‘ k " x av 

tunnel^ are given in figures 8 to 10 . In the case of the sharp- 

edge gu3t, the value of b, chosen for the forcing function Ate 
was determined from the average of values of the time required to 
obtain peak acceleration of the model in the rigid-wing condition 
on entering this type of gust and corresponds to 4.6 chords of 
travel into the gust since, in this problem, the lag in development 
of lift of a wing entering such a gust has the effect of a gust 
with short gradient distance. 


-bt 


Calculations for the model with an arbitrary impressed load 
corresponding to a /ln a of 2g were then made to determine time 

histories of wing-tip acceleration increment, fuselage acceleration 
increment, and wing-tip deflection increment for the conditions 
given in table I. The shapes of the calculated time-history curves 
of fuselage acceleration increment were compared with the experi- 
mental curves for each flight by adjusting the maximum calculated 
value to agree with the maximum experimental value. By using as 
the adjusting factor the ratio of maximum calculated to experi- 
mental values of fuselage acceleration increment for each flight, 
the calculated time histories of wing-tip deflection increment 
were then compared with the experimental time histories. Sample 
comparisons are given in figures 13 to 15 • Included in figures 13 
and 14 is a line representing the calculated value of the static 
deflection 6^ for each condition shown. The calculated wing-tip 

deflection per g of fuselage acceleration increment for each 
condition is compared with the test results in figure 12. 


Additional calculations for the two flexible conditions of the 
model were made under the assumption of no vibratory damping. The 
method used for these calculations was obtained by eliminating the 
damping terms from the left-hand sides of equations (l) and (2) and 
by using net forcing functions which included the damping of the 
vertical motion of the airplane as a whole. The results of these 
calculations are also presented in figure 12 as the wing-tip 
deflection per g of acceleration increment. 


MCA TN No. 1320 


15 


Calculations for Effect of Change of Certain 
Parameters on Dynamic Response 

The airplanes chosen for the calculations are designated 
models A, B, C, and D: model A is a scaled-up version of the test 

model "which may he considered to represent a four-engine landplane, 
models B and C are four-engine landp lanes, and model D is a large 
himotored flying boat. Some of the data for models C and B have 
“been presented previously in reference 11 hut have heen included 
herein for purposes of further analysis. 

The changes in airplane constants chosen for consideration vere 
those resulting from the effects of varying the wing stiffness, of 
arbitrarily omitting all vibratory damping and fuselage clamping 
alone, and of varying the f' Lght conditions of forward velocity and 
weight. The conditions and basic constants for the calculations 
are included in tables II and III. 

General condition s.- The chordwice velocity distribution of a 
single gust va3 assumed to be of the type shown in figure 2(a) . 

In almost all ca.ses the calculations were made for three gust 
shapes, 0-, 10-, and about 20-chord gradient distance, although 
flight experience has shown that, for all sizes of modem airplanes, 
the most probable severe gust has a gradient distance of 10 chords. 

In all ca.ses the forcing function representing the gust was of the 
form Ate'k''. 

Calculations of the effect of repeated gusts were made for 
distributions of the type shorn in figure 2(c) by the method of 
superposition indicated in reference 8, pages 42-1+3, and illustrated 
figure 2 . The values of and H 2 were those chosen for the 

calculations for single gusts and the value of was determined 

to produce the greatest wing-tip deflection for the combination of 
gusts. 

The results of the calculations for each airplane are presented 
in the form of three ratios plotted against the gradient di stance of 
the gust for which the calculations were made (fig. 1 6) . The 
maximum acceleration increments of the fuselage and of the wing tip 
of the airplane in the flexible condition are given as ratios to 
the corresponding maximum acceleration increments when the airplane 
is considered as a rigid body and are called the fuselage acceleration 
ratio and the wing- tip acceleration ratio. The wing-tip deflection 
ratio or the dynamic -stress ratio 5^ mqy /5 s t is the ratio of the 

maximum dynamic wing-tip deflection to the static wing-tip 
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deflection b~ t which would "be determined under static leading 
conditions or as defined previously in equation (3) . The individual 
calculations are outlined in the following sections: 

Model A.- Model A is the scaled-up version cf the gnat-tunnel 
modell The*~two conditions for calculation (tables II and III) 
differ in the same way as the two conditions for the gust-tunnel 
model; that is, the wing frequency for condition 1 is approximately 
half that of condition 2. The results of the calculations in ratio 
form are given in figure l6(a) and 16(b) and a comparison of the 
dynamic-stress ratioe for the two conditions are shown plotted 
against gradient distance in figure 17. 

Model B .- The calculations for model B were made primarily to 
investigate’ the effect of simplifying the calculations by elimi- 
nating the fuselage damping term from equation (2) while beeping 
all other constants identical. The results of the calculations 
for each condition are given in figure 16(c) and 16(d) and a. com- 
parison of the dynamic-stress ratios for the two conditions is 
given in figure 13. 


Model C.- Model C is a modem, large, and fast landplane. 
Calculations were made for a range of gust gradient distances for 
the airplane flying at its cruising speed of 260 miles per hour 
and the results in ratio form are given in figure l6 (e) . Addi- 
tional calculations were made for assumed speeds of 200, 3C0, 
end hOC miles per hour for .he standard gust ■'ri.th gradient distance 
of 10 chords in order to determine the effect of forward velocity 
on the various ratios. The ratios determined, together with those 
for the gust with gradient distance of 10 chords at ?60 miles per 
hour, are shown plotted against forward velocity in figure 19* 

The calculations for the three gradient distances at 260 mile3 per 
hour are used subsequently in obtaining the response to repeated 
gusts . 

Model I).- The calculations for model D were made to show the 
ef f ect~ on" the ratioe of a change in weight cf the airplane such as 
to change the frequency of the fundamental mode of vibration of 
the wing beam.. The two conditions shown in tables II and III for 
this aim lane represent the normal gx’0S3 weight and the overload 
gross weight. Note that the forward velocities for the two con- 
ditions are different. The results of the calculations are given 
in figure l6(f) and 16(g) . A comparison of the dynamic -stress 
ratios for the two conditions is given in figure 20. As in the 
case of model C, these calculations were also used for the deter- 
mination of the response to repeated gusts. 
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Repe ated g usts.- The "basic curves chosen for the extension of 
the calculations to determine the response to two successive gusts 
were the time histories of the reactions of model C, condition 1, 
and model D, conditions 1 and. 2, (figs. 21 to 23) . Sample time- 
history curves for a repeated gust., composed, of two equal and 
opposite gusts of gradient distance of 10 chords arranged relative 
to one another to give maximum negative wing deflection, are given 
in figures 24 to 26. 

Since the interest is primarily in dynamic stress rather than 
in accelerations, the calculated basic curves were superimposed to 
determine the maximum overstress from the combination of the 
reactions to two gusts . The maximum values of dynamic-stress ratio 
occui'red when the repeat gust was negative with respect to the 
fir3t gust and the sequence period or distance (see fig. 2) 

had a pronounced effect on the result. The maximum value of An r 

occurring in the whole sequev.ee was used to determine the static 
wing-tip deflection. Table XV presents selected cases which were 
the most serious of a number of combinations examined. In this 
table, the results refer to repeat gusts having the same velocities 
as the initial gusts. 


DISCUSSION 

Experimental Results and Associated Calculations 


Exp erimenta l results . - It was intended to determine the 
dynamic -stress ratios &dm3x/&st directly from the test results of 
the two flexible and one rigid condition of the model. Examination 
of the test data shoved, however, that the scatter of the data 
within an individual test condition precluded such a procedure. 

As a result, the relative magnitudes of the wing- tip deflection in 
the flexible condition were used for comparison with calculations 
and these magnitudes were determined, as the ratio of the maxiimim 
wing-tip deflection to the maximum fuselage acceleration for a 
given flight of the model. 

As previously pointed out, the experimental data were corrected 
for variation in f onward velocity and gust velocity and., therefore, 
the results for a given gust shape and model condition theoretically 
should be equal. Examination of figures 8 and 9 shows, however, 
that this equality does not exist and that the scatter, sey, in the 
left wing-tip deflections me be a3 great as 17 percent in even these 
selected cases. In addition, the maximum values of acceleration 
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increment for the rigid, condition, from which the static wing-tip 
deflection 6 g .^ would he determined, had large scatter. Since the 
overstress or understress in these cesee would probably not exceed 
20 percent, the experimental data is not sufficiently consistent 
to predict the dynamic-stress ratios • 

Although the experimental data as a whole does not appear to 
he sufficiently accurate to he used as a check of the prediction 
of the dynamic -stress and acceleration ratios hy the analytical 
method, the data in the form given in figure 12 appears to he 
sufficiently consistent to check the order of magnitude of the 
predicted wing-tip deflection. This experimental data in figure 12 
is presented as the ratio of maximum wing-tip deflection to maximum 
measured acceleration increment as a function of the gradient 
distance of the imposed gust . The data for the right wing appeared 
to he less accurate than that for the left wing, probably because 
of local buckling of the right-wing hinge. The test data for the 
left wing then appears to he sufficiently consistent for use as a 
check of the order of magnitude of the predicted wing deflections. 

C ompart son of calculated result s with experimental results . - 
The calculated results for the test model, represented hy solid 
lines, are compared with the experimental results in figure 12. 

The comparison indicates that, for the consistent, experimental 
data for the left wing, experiment and calculation are in good 
agreement as to the magnitude of the deflections. 

Although figure 12 shows that the results for the left wing 
in the model condition of f = 26.1 cycles per second for the 
longest gradient distance differ hy about 15 percent, +his apparent 
disagreement Is minimized hy consideration of the precision of 
measurement of the experimental data . Examination of figure 14(c) 
shows that the recorded wing-tip deflection was only of the order 
of 0.07 inches. Since the precision of measurement is tO.Ol inch, 
the agreement between calculated and experimental results for this 
test condition may he within 5 percent. In general, the average 
of the experimental data for the left wing, with optimum inter- 
pretation of the precision, is lees than 5 percent from the 
calculated values. 

Examination of the time histories on figures 13( a ) and 14(a) 
shows that the actual oscillations of the wing, subsequent to peak 
ving-tir* deflection, a^e not checked hy the calculated reeults. 

This difference is caused hy the necessity of choosing a constant 
division of damping coefficient in the equivalent system to repre- 
sent a given gust condition. The method of application and the 
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amount of vibratory damping, however, appears satisfactory when 
compared in figure 12 with the results obtained when all vibratory 
damping is omitted from the calculations. 

The results of the comparison of experimental and calculated 
data for the flexible wing model then indicate that the method of 
calculation is adequate for prediction of the ratio of maximum 
dynamic wing-tip deflections to maximum fuselage acceleration on 
encountering certain single gust shapes. 


limitations of Calculations 


Limitations shown fry model testa .- The applicability of the 

forcing function Ate“^^ is "best shown by the comparisons of the 
calculated and experimental time histories of acceleration incre- 
ment for the model in the rigid condition (fig, 15) . Since the 
comparisons shown were made by adjusting the calculated curves so 
that the maximum values of acceleration increment agreed with the 
experimental results, any discussion must be based on the shapes 
of the curves alone. It is thought that the greatest part of the 
discrepancies between the curves in figure 15 may be ascribed to 
the effects of pitching motion. The form of the forcing function 
for each gust gradient is intended to include the amount of 


pitching motion indicated by the band of data of shown in 

figure 11. The ratios for the rigid-wing condition of the 

model do not all fall within this band of data. The direction in 
which these ratios differ would tend to explain the difference 
between the experimental and calculated curves of figure 15 , except 
for the long gradient- distance case (fig. 15 (c)) where the experi- 
mental and calculated curves appear to agree perfectly, although 
£0/57.3 

ratio for this case is the furthest from that expected 


the — 


U/V 


(fig. 11) . If, however, the experimental values of acceleration 
increments were raised to consider the pitch correction and the 
calculated values were increased in proportion to the new maximum 
acceleration increment, the two curves would again closely coincide. 
It would appear ^hen that the forcing function in the form 

of Ate is adequate for calculation purposes; however, care 
should be exercised in its application to be sure that in a gust 
the stability of the airplane intended for calculation approximates 
indicated by the band of data in figure 11. Xf not, recourse 
should be made to alternative methods of calculation mentioned 
previously. 
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Ths comparison in figure 


11 of the 


A9/57.; 

u/v 


ratios for the 


flexible conditions of the test model with those for the rigid 
conditions indicates that the assumption of the same pitching 
stability for both rigid and flexible conditions is in error. It 
appears therefore that an investigation should be 3EC.de to deter- 
mine the importance of the effect of wing flexibility on the 
stability of an airplane in a gust. 


Whereas the effect of a constant division of damping coeffi- 
cient appears to be negligible when the present method is used 
to predict the maximum responses to single gusts, the obvious 
overdamping cf the vibrator}- wing motion subsequent to maximum 
ving-tip deflection might well lead to an error, if the results 
for single gusts were superimposed in order to obtain results for a 
succession of gusts. Considerable care should be exercised, there- 
fore, when interpreting results for repeat gusts that are deter- 
mined from single-gust results based on the calculation method of 
this paper. 

In general, the limitations of the method of calculation 
brought out by comparison with results of model teats are thought 
not to affect seriously th 9 results calculated for the response to 
single gusts. A certain amount of dispersion from the pitching 

— *V\4- 

stability assumed by the forcing function Ate does not appear 
to affect seriously the results when presented in ratio form. It 
is apparent that, when the results for single gusts are super- 
imposed to obtain results for successive gusts, the calculations 
for the single gusts must be made with greater attention to factors 
such as the shape of the forcing function and the effect of the 
assumption of a constant division of damping load. 


Oth er llmitat j on s . - Other limitations of the method of calcu- 
lation are apparent In the assumptions that only the fundamental 
mode of wing bending is excited and that the gust is uniform along 
the span of the airplane. When the actual spenwise distributions 
of gust velocity shown in reference 5 are considered, it is apparent 
that tho response in modes other than the fundamental mode of 
bending and the effect of roll and yaw of the airplane should be 
investigated . 

Although, for further simplification, wing torsion was 
neglected in the present analysis, twisting of the wing under 
dynamic conditions in combination with wing bending may readily 
have an adverse effect on the loads for which an airplane must be 
designed. 
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The method of analysis presented in this paper should either 
he extended or replaced by another method so that wing torsion, 
other modes of bending, and the effects of unsymme tri cal gusts 
are considered. 


Results Calculated for Effect of Change of 
Certain Parameters 

Gene ra l resul ts . - The results of the calculations given in 
figure l6 for a.ll the airplanes chosen 6hov that the dynamic-stress 
ratio in all cases increases from an understress for a gust of long 
gradient distance to moderate or high overstres.s for a gust of 
short gradient distance. The general shape of the curves seems to 
indicate that even higher overstress would occur in sharper gusts 
than those gusts examined. The lag in the development of lift in 
the gusts of short gradient distance would preclude, however, such a 
result, since even for on infinitely share gust, the forcing function 
would he similar to that for a gust with a gradient distance of about 
4 chords. Note that, although the dynamic -stress ratios appear to 
approach zero as the gradient distance of the gust increases, this 
condition is not the actual case but results from the fact that, 
when the calculations were r de, account ms not talcen, in the 
division of the damping coeificient, of the changing relative 
significance of the vibration and. of the over-all vertical velocities 
of the wing. If the correct division had been made for the gusts 
with a gradient distance of 2G chords, the dynamic -stress ratio 
would tend to approach 1.0 for most of the cases shown in figure 1 6 . 

The results of the calculations shown in figure 1 6 for models A, 
0, and D shew that, at the design gust of gradient distance of 
10 chords, which is assumed for most conventional airplanes, the 
dynamic -stress ratio varies from about 8-percent understress to 
about 12 -percent overstress. Although at this time a designer 
cannot take advantage of small amounts of understress, when over- 
stress is indicated by the calculations, it is thought that this 
indication should be considered in the design of the airplane. 

Although, for the airplanes considered, the fuselage acceleration 
ratio dees not appear to vary much from a value of 1, the wing-tip 
acceleration ratio at 10 chords is as high as 2.7 in one case and 
greater than 1.6 in most of the cases. This variation indicates 
that the wing-tip acceleration ratio should be examined when con- 
centrated masses or wing components near the wing tip are considered 
in a design. 
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Effect of a cha nge of vine' frequency caused fry a change in wing _ 
stiff ness.- The calculations made for the gust-tunnel tests discussed 
previously 3 how the effect on the dynamic-stress ratios of a change 
of stiffness of a wing such as to halve approximately the frequency 
of vibration of the wing. Examination of figure 17, which gives a 
comparison of the dynamic -stress ratios for the two conditions, 
shows that at present noth:' 3 can "be concluded from the results 
because, although the ratio for the higher wing-frequency case is 
the higher in the shortest gradient gust, the ratio becomes lower 
than the lcw-frequency case as the gradient Increases. At the 
critical gradient distance of 10 chords, the high-frequency case 
shows, however, a reduction in dyne-mi c- stress ratio of 14 percent 
below the low-frequency case. Further analysis of this particular 
Question is therefore needed before the conclusion can be reached 
that a reduction of wing frequency in this manner tends to increase 
the dynamic-stress ratio at the critical gradient distance. 

Effect of simplification by ■ om itting parts of the damning . - 
The comparison of the dynamic-stress ratios for the two calcu- 
lation conditions for model B (fig. 18) and the effect (fig. 12) 
of eliminating vibratory damping from the calculations serve to 
illustrate that simplification of the method by omitting parts of 
the damping does not appear feasible. 

Eff ect of a ch ange in forwa rd velocity of an airplane . - The 
results given in figure 19 illustrate the change brought about in 
the three ratios by increasing the forward velocity of model C 
from 200 to 400 miles per hour. The increase of velocity together 
with the corresponding increase in the rate of application of the 
gust load would appear to result in an increase in the dynamic- 
stress ratio. Figure 19, however, shows that the ratio does not 
vary much as the speed increases and this lack of varis-tion is 
thought to be caused partly by the fact that the aerodynamic 
damping increases directly as the speed and tends to offset the 
expected increase in dynamic- stress. While the results for the 
fuselage acceleration rati show a similar trend, the wing-tip 
acceleration ratio increases from about 1.8 to 2.5 ns the speed 
is doubled, which further emphasizes the recommendation that the 
wing-tip acceleration ratio be considered when concentrated masses 
or wing components near the wing tip are considered in a design. 

E ffec t of a c hange of wing freque n cy caused by a change in 
f light-we ight condition .- The calculations for model D were made 
to show the effect on "the dynamic -stress ratios of a change in 
flight condition from normal gross weight to overload gross weight. 
Table II shows that the forward velocity is different in the two 
cases, but consideration of the foregoing discussion may Justify 
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the assumption that this factor is negligible. The results are 
given in figure 20, with dynamic -stress ratio plotted as a function 
of the gradient dirtance of the gust. A reduction in wing frequency 
brought about by the addition of mass is shown to result in an 
increase in the dynamic -stress ratio. Further analysis and calcu- 
lations appear to be needed, but a tentative conclusion is made 
that for a conventional airplane a reduction in wing frequency 
caused by the addition of mass results in an increase in the 
dynamic -stress ratio. 


Effect of repeated g u r ' 3 on the dynamic -stress ratios . - The 
results of calculations fc.u models C and D were presented previously 
in reference 11, and the values of table TV are talren from table IX 
of that paper and show the most serious values of overstress for 
two gusts from a combination of the calculations for a single gust. 
Table IV shows that no definite correlation exists between the effect 
of gradient distance of the first end second gusts and the distance 
between them This lack of correlation results from the 

influence of certain other factors, such as the relation between 
the time to peak acceleration and the period of wing vibration, 
which complicate the problem when the reactions to one gust are 
superposed on the reactions to another gust. 

Examination of the values in table TV shows that substantial 
overstress exists for all the combinations of gusts presented and 
that the addition of a short gradient gust produces the largest 
valxie of overstress. As indicated previously, however, the gust 
velocity measured by an airplane tends to decrease from a maximum 
for a gradient distance of 10 chords os the gradient distance is 
decreased. As a result, the values of dynamic stress to be con- 
sidered are those (indicated by the footnote in table IV) that 
represent the combination of the reactions to gusts heving a 
gradient distance of 10 chords. Before an estimate can be made 
as to whether the overstress shown for the airplanes in question 
is serious, it is necessary to consider the conditions upon which 
the values in table IV were calculated and the effect of the 
intensity and size of gusts end their spatial distribution in the 
atmosphere . 

The results shown in J able IV were based simply on the premise 
that the quantity of interest was H he ratio of maximum stress 
obtained under dynamic conditions to the maximum stress that would 
be computed under static conditions without regard to velocity or 
spacing of such pairs of gusts in the atmosphere. Examination of 
the time histories of res-ctions given in figures 2k to 2 6 indicates 
that for each case shown the dynamic -stress ratio given in table IV 
is the ratio of the maximum dynamic deflection in the second gust 


2b 


MCA TN No . 1320 


to the static deflection 6 gt computed from the maximum value 
of Axij. in the second gust. As pointed out in refei’ence 11, 
however, when design conditions are considered, the ratio to use 
is that of the maximum dynamic deflection of the whole sequence 
to the deflection computed from the static load in the first gust 
"because, if the two grata w :a each of design velocity, the 
design static load would "be attained on passage through the first 
gust. The values determined, on this basis for the succession of 
two gusts of 10 -chord gradient distance indicated in the table 
would then be as follows: 


Model 

Condition 

5 d/' 5 st 

C 

1 

1.62 

D 

1 

I.58 

D 

2 

1,6b 


An analysis given in reference 11, however, based on frequency 
data of single gusts, indicated that the gist velocities of two 
repeated gusts would range from 0.6l to 0,75 the velocity of the 
single design gust so that, multiplied by these ratios, the 
dynamic-stress ratios given would be reduced to an average 
of 1.10 times the design stress. 

A recent statistical analysis of the characteristics of 
repeated gusts in turbulent air (reference 12) provides more 
concrete data, however, than the analysis used in reference 11. 

Two conclusions from reference 12 state that sets of two repeated, 
gusts with average absolute effective gust velocities of 25 feet 
per second apparently are encountered in turbulent air as often 
as single gusts of intensity greater than 30 feet per second, and 
that the over-all average c .icing between two repeated gusts is 
about 25 chords. Note that spacing as defined in the reference 
paper i3 the distance between acceleration peaks and in the 
terminology of this paper would correspond to the sum Eq + H^. 

Table TV shows that the sums of Hg + in the cases 

indicated by the footnote do not approximate 25 chords; therefore, 
a new superposition of the responses to single gusts of 10 chords 
was made so that the suacing would be 25 chords. The dynamic- 
stress ratios for the three cases were then determined as the 
ratios of the maximum dynamic deflection in the sequence to the 
deflection computed from the static load imposed by the first gust 
and the results were reduced by the ratio of the average gust 
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velocity for two repeated gusts to the design velocity for a single 
gust which is 25/3 0. The results determined on this "basis are 
compered with the results for the single gusts of 10 chords in the 
following table : 


Model 

Condition 

Design 

repeated gusts 

Design 
single gusts 

5 W /6ct 

C 

1 

0.96 

1.07 

D 

1 

1.08 

0.92 

D 

2 

1.10 

1.09 


The variation in the results indicates that the dynamic-stress 
ratios for a design repeat gust should "be investigated. 'The 
results also indicate, however, that these dynamic-stress ratios 
are not likely to be much greater than those which would be 
determined for a design single gust. 


CONCLUSIONS 

Analytical Method and Experimental Work 


It appears from consideration of the comparison of the experi- 
mental work and associated calculations that: 

1. The analytical method as presented in this paper is of 
sufficient accuracy to predict the ratio of the maximum dynamic 
wing-tip deflection increment to the maximum fuselage acceleration 
increment for a conventional airplane. 

2 . Simplification of the method by omitting parts of the 
damping does not appear feasible. 


Calculations 


The analysis of the results of the calculations for the effect 
of change of certain gust and airplane parameters indicates : 

1. Tne dynamic -stress ratio for airplane wings from encountering 
gusts increases as gradient distance decreases. 
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2. For the assumed design gust with the gradient distance 
of 10 chords, the overstress in a single gust may he as much 
as 12 percent. 

3 . Although the results for the two airplanes each of which 
encountered one type of repeated gust were not considered con- 
clusive, a repeated gust does net seem to he mere critical than a 
eingle gust. 

H. For the assumed design gust of 10-chord gradient distance, 
the wing-tip acceleration increment is generally much greater than 
the fuselage acceleration Increment and should he taken into 
account when designing for concentrated masses or wing components 
near the wing tip. 

6. For the assumed design gust of 10-chord gradient d.i stance 
the dynamic -stress ratio do< not change appreciably with change 
in forward velocity oi ‘.he , irplane, hut an increase in speed is 
accompanied hy an increase in the wing-tip acceleration, ratio. 

6. A reduction in wing frcauency either hy a change in stiff- 
ness or hy a change in weight increases the dynamic -stress ratio 
at the assumed design gust of 10-chord gradient distance. 


Langley Memorial Aeronautical laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., January 15, 19I+7 
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APPENDIX 


DETAILS OE METHOD A Iff) PROCEDURE FOR 


DETERMINATION OF CONSTANTS 


Certain as Gumptions were made to simplify the problem of 
determining the dynamic response of airplane wings on encountering 
gusts. These assumptions are: 

(1) The imposed gust loads are symmetrical about the center 
line of the airplane and their characteristics in the line of 
flight are known . 

(2) The imposed gust loads excite only the fundamental mode 
of bending of the wing with the result that the stress in the wing 
is proportional to the deflection of the wing tip. 

( 3 ) The forward velocity of the airplane is constant during 
passage through the gust. 

Together with these assumptions, an airplane is reduced to 
the equivalent biplane shorn in figure 1. The equations of 
motion of the two parts of the equivalent biplane are reduced to 
linear equations with constant coefficients. The solution of the 
equations and the method of determination of coefficients follow. 


General equations .- Tho equations for which the solutions 
follow contain the type of forcing function used for calculations 
in this paper. 


Solutions of Equations 



(Al) 







K(5v - 6f) = Af 0 te 


(A2) 
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Since equations (Al) and (A2) are simultaneous linear 
differential equations, equation (A2) may "be solved for 8 W and 
the result substituted for b v in equation (Al) . This process 
leads to the following: 


(V + a-jl)3 + a.pD^ + a 3 D)8 f = P 1 te' 1>t + F 2 * 


-ht 


(A3) 


where 


X, t Mf> -4- M, r 
•‘•e ± e ~ ' T e 

e l “ M M- 

v e f e 


(A.Ua) 


a o = 


K(Mw e + MfJ + X Wf Xf 


e e 


M M f 
I e 


(A 4b) 


e 3 


+ *0 


(A4c) 


F, = 


A Wp K Af 


1 “ M v Mf Mf 

© 0 G 


7e J 


Mw e T M v 


(A5a) 




(A5b) 
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Equation (A3) may "be recognized as a linear differential 
equation with constant coer . J.cients which may he solved hy methods 
commonly used for solving differential equations. Solution 
for is 


Of = 


-Rit 


c l e 


(cos R 2 t) 


-R-it . . -Rot 

c 2 e x (sin Rgt) + c-e -> 


+ + 



+ 



(A6) 


where R-^ and P.p are the real and imaginary parts of the complex 
roots P]_ + Rg\/-1 and R, is the real root of the cubic equation 
derived from equation (A 3 ); c-j, c 2 , c-^, and are constants of 
integration; and. in the particular integral of the solution; 

G, = b^ - a,b3 + a 0 t>2 - a,b (A7a) 

x id 3 

G £ = Vb3 - 3 &1 h 2 + 2a n b - a^ (ATb) 


In the case- where G-^ = 0 a special solution of the equation 
would he required, hut this solution can he obviated hy choosing a 
slightly different value for h. 

Equation (A?) is solved for the deflection of the equivalent 
wing with respect to the equivalent wing-fuselage (5 v - 6^) and 

for the space position of the equivalent wing 5 W . The resulting 
equations follow: 


6 


w 


Sf 


-££ - Af e te " V 

~K dt 2 + K dt K 


(A8) 


w 
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A fe te 
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K 


(A9) 
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The vertical velocities of the equivalent wing-fuselage and 
of the equivalent wing are determined from the first derivatives 
of equations (a 6) and (A9) , respectively, and the normal accelera- 
tions are determined from the second derivatives. 

In order to perform an arithmetical solution for a particular 
airplane, it has been found advisable to insert the pertinent 
numerical constants for the airplane in equation (A3) and. 3 olve . 

In addition, experience has shown it necessary to carry eight to 
ten significant figures throughout the solution so that the 
results are useful "because, at different stages in the evaluation, 
small differences of large quantities are obtained. 

Rigid -wing acceleration increment .- The normal gust -load 
design procedure assumes that the load is applied statically and 
that the normal acceleration is constant along the span of the 
airplane. The difference "between this assumption end the actual 
case is shown by ratios of the accelerations determined under 
dynamic conditions to the accelerations that would be determined 
if the airplane were rigid. In order that the comparison be valid, 
the rigid airplane acceleration Increment An r is determined as 

the response of the equivalent biplane to the over-all forcing 
function on the airplane when the springs are replaced by rigid 
rods. With this restriction equations (Al) and (A2) may be 
combined to become 


0 t 

Aiij. + \ c i An r 

Jo 


dt = A c te"^ 


(A10) 


where 


and 


a 2 s w d 2 $f 


An « — 1 r 


dt 


dt ' 
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When the integral in equation (A10) ia removed "by differ- 
entiation, the equation he come a a linear differential equation 
with constant coefficients which may he solved in the same manner 
as equation (A 3 ) . The complete solution is: 


An r = 


A. 


\ r - h 


r- + 


A c h 


( X c - 


W 




- e 


-X,t 


A-r 


A c l. 


te 




(All) 


S tatic deflectio n end dynani c~ stress rat i o . - The static 
deflection used in the dynamic -stress ratio to represent the 
static stress in the wing is considered as the deflection of the 
upper wing of the equivalent biplane under the conditions of the 
normal static design procedure. The static deflection is then 
computed as follows: 


b et 


_ ^ r max( LiW ' Ww e) 


K 


(A12) 


The ratio of the raximum value of dynamic wing-tip 

deflection 5, , as determined hy equation (A8), to the static 

max 

deflection 5^ is called the dynamic -stress ratio since the 
stress, determined under static conditions, is multiplied by this 
ratio to take into account dynamic conditions. 


Determination of Constants for the Equations 

Req uired condition s . - In determining the values of the coef- 
ficients for the equivalent biplane, the following conditions 
should he satisfied: 

(1) The total mass and the total load on the equivalent biplane 
should he identical with those of the original airplane. 

(2) The kinetic energy of vibration of the upper wing should 
closely approximate that of the original wing beam for an amplitude 
of vibration of the upper wing eaual to that of the tip of the 
original wing. 

(3) The natural frequency should be the same as that for the 
fundamental mode of the original wing. 
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(4) The upper wing should deflect under the equivalent static 
load the sane amount that the wing tip of the airplane does under 
the corresponding aerodynamic static-lead distribution. 

(5) The damping coefficient of the upper wing should represent, 
at least to peak load, the damping of the motion of the wing. 

Equivalent masses and Mp .- The equivalent mass M-. 

is obtained from an approximate requirement that the kinetic energy 
of vibration of the original wing beam is reproduced by the upper 
wing of the equivalent biplane. The equivalent mass of the 
fuselage Mf e is taken equal to the total mass of the airplane 

minus M w . 

w e 

In the absence of mor definite information the combined 
effect of concentrated manses and nomuiiform wing structure is 
assumed to be such that, the equivalent wing mass may be determined 
from the following relation which represents an approximation to 
that which would be derived for a uniform cantilever beam. Thus, 





(see fig. 1) where x is the distance from wing root to tip and x 

is the distance to the individual mass AM. For special cases, such 

ps strut-supported wings and arrangements other than the conventional 

cantilever monoplane, other suitable approximations can be devised. 

If the wing were a uniform cantilever beam, ML would be approxi- 

e 

mately 2h percent of the mass of the wing (reference 13, PP* 83 - 89 ) • 

Equ ivalent spring constant K .- The equivalent spring constant 
defines the springs in the equivalent biplane which allow approxi- 
mately the same frequency as the cantilever wing of the original 
airplane. The value of K may he approximated by using M w and 

the known or estimated wing frequency in the relation: 
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Air loa ds or forc in g func tj.or g.- As in the case of the mass of 
the airplane, the total air load must he divided into the components 
affecting the motions of the fuselage and wing tip of the original 
airplane. If the deflection of the wing tip under a given static 
air load is known, the transfer to the equivalent "biplane is 
accomplished ty merely dete ..lining the load applied to the equivalent 
spring which gives that same deflection to the upper wing of the 
equivalent biplane. The remainder of the load i 9 then applied to 
the lover wing-fuselage combination. Any air load of the same type 
as this static air load may be divided into the same proportion and 
applied to the equivalent system to give the same deflections and 
reactions as the load causes on the original airplane. 

Data on deflection, however, is often not available. Recourse 
may then be made to several methods of determining the division of 
loads. One simple approximation is to assume that the wing of the 
airplane is a uniformly stiff cantilever from root to tip. Then, 
if the beam is considered to be weightless, the equation for the 
deflection of the tip of the beam under a concentrated load placed 
at the tip is compared with the equations of deflection for two types 
of loading distributed along the beam as follows: a uniform load 

along t he beam and a loan uniformly tapering from a maximum at the 
root to zero at the tip. The differences in the three equations 
ere in the numerical "efficiency" factors which are one-third for 
the concentrated load, one-eighth for the uniform load, end one- 
fifteenth for the tapered load. A concentrated load at the tip 
of. a cantilever beam that will give the same deflection as a uni- 
form lead is then one-eighth divided by one-third or three-eighths 
of the un.ij.orm load, and the equivalent concentrated load to replace 
the tapered load is three -fifteenths of the tapered load. The esti- 
mate Cj. the shape of the span leading of a given airplane probably 
falls somewhere between that of a uniform load and of a tapered 
load so that a comparison of these shapes indicates that between 20 
and 37.5 percent of the total load produces the tip deflection and 
the remainder is considered as acting at the wing root. If the 
wing under consideration does not approximate closely enough a 
unnorm cantilever, allowance for this discrepancy can be made in 
choosing the percentages of load. 

The forcing functions A„ te - ' bt and A* te“‘ bt in 

w e x e 

equations (Al) and (A2) represent the air load on the wing apart 
from the damping, divided as in the previous discussion, where 




= A 
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and 
A = 


An*, W 

at lit = 1 or WbeAr^ 


(A13) 


te 


-bt 


An arbitrary load factor Increment the airplane wu Id experience 
a if it had no vertical motion when traversing the gust 


b = - at bt = 1 or at the maximum value of the function te 
t 


-bt 


The next problem is the determination of a value of b so 
that the forcing function represents a gust of a given shape. 
Gust-tunnel tests of relatively stable end rigid models { reference 6 ) 
have shown that the peak of the acceleration- increment curve occurs 
at the end of Hj_, as given in figure 2(a). The time necessary 

to reach the maximum value of the time history curve of the rigid- 
wing acceleration increment, which is obtained from the forcing 
function curve by use of equation (All), is assumed, therefore, to 
represent the length of the gradient of the gust imposed on olie 
airplane. Since the peaks of the forcing -function curve and the 
acceleration-increment curve do not occur at the same time except 
in the case of an eirplane of infinite mass, the following pro- 
cedure for the determination of a value for b has been derived. 


(1) Three or four va” es of b are chosen so that they 
represent a range of peak values of the forcing function corre- 
sponding to from 2 chords of travel of the airplane to 40 or 50 chords 

of travel. 

(2) The time necesaary to reach the maximum value of An r is 
determined for each value of b chosen and converted to chords of 
travel of the airplane . 

( 3 ) A plot is made of b against chords of travel to peak 
of the An r curve. 

(4) The values of gust gradient distance chosen for the calcu- 
lations are then used with the curve plotted in step 3 to determine 
the corresponding values of b. 

If the forcing function Ate _tt does not adequately represent 
the case considered, the actual forcing function for the sharp-edge 
gust may be determined from calculations similar to those calcu- 
lations outlined in reference 7 or by recourse to gust-tunnel tests. 


MCA TN No. 1320 


35 


Equations (Al) and (A2) would then "be solved for a unit -function 
type of forcing function (> Terence 8 ) and the result built up into 
the response of the fa.e 3 r.it . 3 airplane to a sharp-edge gust by the 
graphical method outlined in reference 9 • This response of the 
flexible airplane to a sharp-edge gust may then be considered as 
the response to a unit -function -type gust and be again built up to 
represent the response of the flexible airplane to any type of 
symmetrical gust. 

Damping f ac tor Consideration of the reaction of an air- 
plane to a given gust or the reaction of the equivalent biplane 
to the forcing function indicates that two distinct motions have 
to be damped - the vertical motion of the airplane as a whole and 
the vibratory motion of the wing itself. Since the vertical motion 
rises from the action of unsteady lift, the damping of this motion 
is a.lso subject to unsteady lift effects. The vibratory motion of 
the wing falls in this category because the lift of an oscillating 
airfoil has been shewn (references 14 and is) to be affected by 
unsteady lift phenomena. Determination of the effects of this 
unsteady lift is therefore necessary in order to predict correctly 
the dynamic stress of an airplane wing upon entry into a gust. 

The obvious solution would be to include the equations of unsteady- 
lift damping directly in the dynamic -stress equations, but this 
procedure would destroy the linearity of the equations and make 
them very difficult to solve. The steady-lift damping would 

be £ mSY times the velocity of oscillation with £ mSV a constant 
2 2 

for a given cese. Inasmuch es having the damping coefficient in 
this form would fulfill the conditions for normal solution of the 
original equations, an ana!!' .is and some tests were made to 
determine whether the effc . . of unsteady lift on the damping force 
could be considered to reduce the steady-!’ i+’t damning force by a 
constant factor - a damping-efficiency factor - without seriously 
impairing the reeults. 

The analytical determination of this factor was derived frem 
the following expression for the change in Cp brought about by a 
sudden change in anglo of attack: 



(AIM 
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where 

AC-j^ change in coefficient of lift 

m elope of lift curve, radians 

C-^ ratio of absolute value cf the unstesdy-lift function for a 

a sudden change in angle of attack to its absolute value 

at t = 00 for a wing of aspect ratio of 6 

s^' distance of travel from initiation of angle-of -attack change 
to point at which lift is desired, half chords 

s' distance of travel from initiation of angle -of -attack change, 

half chords 


SSL rate of change of angle of attack with half chord of travel 
da' 

Substitution of Jones' approximation of the unsteady -lift function 
for aspect ratio 6 for a sudden change in angle of attack (reference 1 6 , 
equation ( 29 )) changed to ratio form 


is made in equation 
motion of the wing, 


s') = 1.000 - 0.36le‘• 0 •3 8l ( s l , - s, ) 

(Aik) . Then, assuming a simple harmonic vertical 
the following substitution is made: 


where 

U 

a = — 
Y 

U 

Y 

X 


da 

d3 


Xk 


Y - -y- sin ks' 


vertical velocity of wing 
horizontal velocity of wring 

measure of the maximum vertical velocity of wing 


k = 


2rt 
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and 

0 ' distance cf travel of wing for one period of oscillation, 
half chords 

The indicated integration in equation (Al4) is performed and 
the result is divided "by the assx^med steady -lift damping such that 
the damping-efficiency factor is given by 


-C^ cos ks^* - Cg sin ks*j 1 ■+■ 1.000 + C-^e ^ o i 

1 - cos k8]_' 


(A15) 


where 


C-l = l.OCO 


0 . 361k 2 
0.145 + k 2 


0.136k 
0.145 + k 2 


0.361k 2 
, ' 2 
0.145 + k 


Evaluation of expression (A15) for a series of ratios of wing 
frequency to half chords of travel shows that the darcping-eff Jciency 
factor varies considerably throughout a cycle of oscillation of the 
wing, but the greatest variation occurs as the vibration velocity 
approaches zero or at the point of minimum damping force. As a 
result, an average value of 75 percent for the magnitude of the 
damping-efficiency factor was taken from the parts of the cycle 
where the damping force was near the maximum value so that the 
total damping coefficient X is given by the equation: 


X = 0.75 | 
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In order to obtain an experimental verification of the method 
and result of the theoretical determination of damping-efficiency 
factor, tests were conducted on a 10 -inch half -span half -mod el at 


three different wind speeds in the scale model of the full-scale 


tunnel. Oscillatory motion about a flexure plate hinge at the wing 
root was begun by the decrease in lift caused by the sudden deflec- 
tion to zero of a plain flap that formed the entire trailing edge 
of the wing. The flap motion and the wing-tip motion were recorded 
as time histories and the damping coefficients were then determined 
from the logarithmic decrements of the wing-tip oscillations 
(reference 13 , p. 35) • The results of these tests indicated that 
the value of 75 percent determined by the theoretical analysis was 
a valid quantity for the damping-efficiency factor. 


The damping coefficients for the upper wing and lower wing- 
fuselage must now be determined. The proper values for the 
individual coefficients are dependent upon the spanwise distri- 
bution of the vertical velocities of the original wing which 
changes as the airplane penetrates the gust. In the case of 
uniform spanwise distribution of vertical velocity (case 1) , the 
airplane damping coefficient X is divided in the same way as 
the impressed air-load coefficient so that 




X 


For the case where the spanwise distribute on of vertical velocity 
equals the deflection curve of the assumed uniform cantilever wing 
of the original airplane (case 2), the total damping load on the 
original wing is concentrated near the wing tip and is approxi- 


mately equal to X times one -third the tip velocity. When this 
situation is applied to the equivalent biplane, the total damping 
on the equivalent biplane is applied to the upper wing so that X w ^ 

is then equal to -^X . 


Case 2 represents the situation as the airplane first pene- 
trates the gust and case 1 represents the situation later in the 
penetration when the vibration has damped out . In order to avoid 
destroying the linearity of equations (Al) and (A2) , coefficients 
chosen on the basis of conditions existing at maximum load ere used 
to calculate the entire time history. If case 1 exists, the time 
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histories of reactions calculated are not seriously affected^ 
however, when case 2 exists, the time histories subsequent to peak 
load would "be quite seriously in error. In order to account for v 
this inaccuracy and to effect a transition "between the two cases, 

Xf^ is taken equal to -X in case 2. This assumption will not 

affect the results for the penetration into the gust represented 
"by case 2, "because the vertical velocity of the fuselage of the 
original airplane and the wing-fuselage of the equivalent biplane 
would be zero. 

In order to provide a way to determine which division of 
damping coefficient is applicable in a given calculation, an 
arbitrary criterion was found that depended on the relation 
between the airplane vertical velocity, which wa 3 determined by 
integrating the time history of the rigid-wing acceleration incre- 
ment £n r to its maximum value, and the vibration velocity of the 

wing, as determined by dividing the static deflection 6 0t by the 
time required to reach maximum An r from the beginning of the gust. 
When the ratio of airplane vertical velocity to the vibration 
velocity is less than 5, the division into one-third and two-thirds 
applies, and when the ratio is greater than the damping coeffi- 
cient is divided into one-fourth and three -fourths . Note that the 
divisions of damping coefficient given apply in particular to con- 
ventional cantilever-wing airplanes. Other types of airplanes 
would require special analysis. 


40 


NACA TN So. 1320 


REFERENCES 


1. Kussner, Hans Georg: Stresses Produced in Airplane Wings "by 

Gusts. HACA 1M No. 654 , 1932. 

2. Bryant, L. W., and Jones, I. M. W.: Stressing of Aeroplane 

Wings due to Symmetrical Gusts. R. & K. No. l690, 

British A.E.C., 1936. 

3. Williams, D., and Hanson, J.: Gust Loads on Tails and Wings. 

E. &. M. No. 1823, British A.B.C., 1937* 

4. Sears, William B., and Sparks, Brian 0*: On the Eeaction of 

an Elastic Wing to Vertical Gusts. Jour. Aero. Sci., 
vol. 9, no. 2, Dec. 1941, pp. 64-67. 

5. Moskovitz, A. I.: XC-39 Gust Eesearch Project - Preliminary 

Analysis of the Lateral Distribution of Gust Velocity along 
the Span of an Airplano. NACA EB, March 1943 • 

6. Donely, Philip: An Experimental Investigation of the Normal 

Acceleration of an Airplane Model in a Gust. NACA TN Ho. 706, 
1939. 

7. Donely, Philip, Pierce, Harold B., and Pepoon, Philip W.t 

Measurements and Analysis of the Motion of a Canard Airplane 
Model in Gusts. NACA TN No. 758, 1940. 

8. Berg, Ernst Julius : Heaviside's Operational Calculus. 

McGraw-Hill Book Co., Inc., 1936 » 

9. Jones, Robert T.: Calculation of the Motion of an Airplane 

under the Influence of Irregular Disturbances. Jour. Aero. 
Sci., vol. 3, no. 12, Oct. 1936, pp. 419-425 • 

10. Rhode, Richard V.: Gust Loads on Airplanes. SAE Jour., 

vol. 40, no. 3, March 1937, pp* 81-88. 

11. Pierce, Harold B.: Dynamic Stress Calculations for Two 

Airplanes in Various Gusts. NACA ARE, Sept. 1943.. 

12. Moskovitz, A. I., and Peiser, A. M.: Statistical Analysis of 

the Characteristics of Repeated Gusts in Turbulent Air. 

NACA ARE No. L5H30, 1945- 


NACA TN No. 1320 


41 


13. Timoshenko, S.: Yihration Problems in Engineering. Second ed., 

D. Van Nostrand Co., Inc., 1937. 

14. Kiissner, H. G.: Zusammenfassender Bericht uber den instationaren 

Auftrieb von Flugeln. Luftfahrtforschung, Bd. 13 , Nr. 12 . 

Dec. 20, 1936, pp. 410-424. 

1?. Theod-orsen, Theodore: General Theory of Aerodynamic Insta- 

bility and the Mechanism of Flutter. NACA Rep. No. 496, 

1935- 

16. Jones, Robert T.: The Unsteady Lift of a Wing of Finite 

Aspect Ratio. NACA Rep. No. 68l, 1940. 


TABLE I 


•\ 


CHARACTERISTICS OF THE TEST MODEL AND 
EQUIVALENT CONSTANTS FOR CALCULATION 


Values constant for tests 

Weight, lb 

Wing area, sq ft 

Mean geometric chord, ft 

Span, ft 

Slope of lift curve, per radian 
Forward velocity, fps 
Oust velocity, fps 

1. 832 
1.183 
0.394 
?.o 

>73 

61.0 

6.0 

Pitching moment of inertia, slug-ft^ 
Radius of gyration of wing, ft 
Weight of wing, lb 
Mass ratio, % /m 

Air load ratio, A w /M 

Damping factor, X, lb-sec/ft 

0.00782 

O .663 

0.293 

O.O 333 

0.37 

0.2954 

Values changed in tests 


Variables 

f = 13.5 ops; K = 13.608 lb/ft 

f = 26.1 cps; K * 51.024 lb/ft 

Gust 1 

Gust 2 

Gust 3 

Gust 1 

Gust 2 

Gust 5 

Gradient distanoes of 
gust, chords 
Damping ratio, \ w /K 
w e 

Forcing-funotion factor, 
A, lb/sec 
Time constant, b, 
per sec 

4.6 

0.50 

275.948 

27.70 

7.8 

0.50 

151.921 

15.25 

21.2 

0.37 

37.2580 

3-74 

4.6 

O. 5 O 

275*948 

27.70 

8.8 

O .50 

152.495 

13.30 

22. 4 
0.37 

38.7523 

3.89 
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TABLE II 

CHARACTERISTICS OF AIRPLANES CHOSEN 
FOR CALCULATIONS 


Airplane 

Condi- 

tion 

Weight 

(lb) 

Wing 
area 
( sq ft) 

Wing 
loading 
(lb/sq ft) 

Span 

(ft) 

Mean geo- 
metric 
chord 
(ft) 

Natural 
wing fre- 
quency 
( cps) 

Slope of 
lift curve 
(per 
radian) 

Number 

of 

engines 

Forward 

veloc- 

ity 

(mph) 

Remarks 

Model A 

1 

100 , 000 

1700 

58.8 

113.8 

14.93 

2.20 

4.73 

4 

256 

Hypothetical 

airplane 


2 

100,000 

1700 

58. 8 

113.8 

14.93 

4.25 

*-73 

4 

25b 

scaled from 
flexible 












wing model 

Model B 

1 

44,860 

2480 

18.1 

150.0 

18. 53 

3-89 

4.76 

4 

200 

Fuselage 

damping 


2 

44,860 

2480 

18.1 

150.0 

18.53 

3-89 

4.76 

4 

200 

omitted in 
condition 2 

Model C 

1 

100 , 000 

1710 

58.5 

140.0 

12.21 

2.45 

5.04 

4 

260 

Speeds 


2 

100,000 

1710 

58.5 

140.0 

12.21 

2.45 

5.04 

4 

200 

assumed 

for 


3 

100,000 

1710 

5 * 5-5 

140.0 

12.21 

2.45 

5.04 

4 

300 

calculation 


4 

100 , 000 

1710 

58.5 

140.0 

12.21 

2.45 

5.04 

4 

400 


Model D 

1 

62,500 

1826 

34.2 

140.0 

13.04 

2.50 

4.93 

2 

190 

Normal gross 
weight 


2 

102 , 000 

1826 

56.0 

140.0 

13.04 

1.43 

4.93 

2 

160 

. 

Overload 
gross weight 
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TABLE III 


CHARACTERISTICS AND CALCULATION DATA FOR BIPLANE EQUIVALENTS 
OF AIRPLANES IN TABLE II 


Biplane 

equi- 

valent 

of 

Condi- 

tion 

Gradient 

distance, 

( chords ) 

Weight 

(lb) 

M v e 

( slugs ) 

Mw 

M 

Damping 

factor, 

A 

( lb-sec/ft) 

X w e 

or 

Forcing 
function 
factor, A 
( lb/sec) 

-IS 

A 

Spring 

constant, 

K 

(lb/ft) 

Time 

constant, 

b 

(per sec) 

Model A 

1 

4.6 

9.0 

22.0 

100,000 

103.33 

0.0333 

2,614.0 

0.500 

2.445.000 

1.348.000 
329,000 

0.37 

19,630 

4.R0 

2.4g 

.606 

2 

4.6 

g.g 

22.4 

100,000 

103.33 

.0333 

2, 614.0 

.500 

2.445.000 

1.169.000 
342,500 

• 37 

73,300 

4. 50 
2.15 
.630 

Model B 

1 

2.8 

**■•7 

5-7 

7.1 

44, 860 

65.43 

.0470 

3,090.0 

• 333 

975,544 

487,772 

324,368 

243,886 

.25 

56,436 

4.00 

2.00 
1-33 

1.00 

2 

3-3 

6.7 

8.5 

10.0 

44,860 

65.43 

.0470 

1,030.0 

1.000 

975,544 
487,772 
324, 368 
243,886 

.25 

56,436 

4.00 

2.00 
1-33 

1.00 

Model C 

1 

3-7 

10.0 

20.0 

100,000 

106.38 

.0343 

2,972-9 

•333 

3,773,000 

1,255,855 

482,226 

•25 

25,233 

6.94 

2.31 

.887 

2 

10.0 

100,000 

106.38 

.0343 

2,256.3 

.333 

989,461 

•25 

25,233 

• 1.82 

3 

10.0 

100,000 

106 . 38 

•0343. 

3,384.5 

.333 

1,440,699 

.25 

25,233 

2.65 

4 

10.0 

100,000 

106.38 

.0343 

4,512.7 

.333 

1,980,114 

• 25 

25,233 

3.64 

Model D 

1 

4.3 

10.3 

20.5 

62 , 500 

50.23 

.0258 

2,239.0 

• 333 

1,450,893 

428,132 

1.55,962 

.25 

12,406 

4.27 

1.26 

.459 

2 

3-9 

10.1 

20.2 

102 , 000 

154.48 

.0487 

1,885.4 

• 333 

2,168,225 

720,893 

273,939 

.25 

12,406 

3-91 

1.30 

.494 
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TABLE IV 

MOST SERIOUS VALUES OP OVERSTRESS PROM ADDING 
THE REACTIONS OP TWO GUSTS 

[Based on reference ll] 


H 1 

(chords) 

h 2 

( chords ) 

H 3 

(chords) 

&d 

a w 

Set 

Model C, condition 1 

3.75 

3-75 

16. 2S 

1.40 

3.75 

9-99 

13.74 

1.26 

3-75 

19- 9S 

3-lb 

1.26 

9*99 

3-75 

39 ..94 

1-35 

9-99 

9-99 

37-46 

a 1.25 

9.99 

19.9S 

27- *7 

1.0s 

19 . 9 s 

3-75 

72.43 

1.44 

19.9S 

9-99 

b9-93 

1.34 

19.9S 

19.9S 

59.94 

1.20 


Model D, condition 1 

• 

4.26 

4.26 

14.50 

1.20 

4.26 

10.25 

11.96 

1.09 

4.26 

19-53 

5.13 

1.09 

10.25 

4.26 

35.04 

1.11 

10.25 

10.25 

32.4S 

a i. 21 

10.25 

19.53 

25-64 

1.20 

19-53 

4.26 

69-22 

1.42 

19.53 

10.25 

66.66 

1.34 

19-53 

19-53 

59-82 

1.39 

Model D, condition 2 

3-96 

3.9^ 

15.47 

1.47 

3-96 

10. OS 

13.31 

1-35 

3-96 

20.15 

11. 87 

1-35 

10. OS 

3.9b 

35-9S 

a 1 * 4 ! 

10. OS 

10. OS 

33. S3 

a 1. 26 

10. OS 

20.15 

32.39 

1.11 

20.15 

3-96 

73-41 

I.56 

20.15 

10. OS 

71.25 

1.36 

20.15 

20.15 

69. SI 

l.?3 


a Used In extended analysis of repeated gusts. 
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Fig ure I.- E g ui valent biplane. 
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Fig. 2 NACA TN No. 1320 


NACA TN No. 1320 


Fig. 3 



Figure 3.- Test model, 
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Figs. 5,6 
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Figure 5.- Diagram of knife edge, strata, and fuselage 
Spring of test model. 




Figure 7.- Instrument record showing A, wing -tip -deflection time 
histories, B, accelerometer time history, and C, distance 
between A and B resulting from the record light beams 
striking the film 90° apart on the drum. 
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U/U max 0 y and fuselage Left wing-tip Right wing- tip 

acceieratfon increment, An ,g c/ef/ection increment, in . defection increment, in. 
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Fig. 8a 



U/f mat a v 




Horizontai d /stance from front edge of gust tunnel , 

chords (referred to eg. of model) NATIONAL ADVISORY 

COMMITTEE FOR AERONAUTICS 

(a) Sharp -ec/ge gust with Z. 4 -chord gradient distance. 

Figure 8 .-History of events . f =/3.ffcps. 
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(b) Gust with 7.8 -chord gradient distance. 


Figure 6. - Continued. 
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(a) Sharp- edge gust with /.8-chord gradient distance. 

Figure 3. -History of events > f- F6. / cps. 
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(b) Gust with Q. 8- chord gradient distance. 


Figure 9. - Continued. 
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(c) Gust with ZL4- chord gradient distance. 

Figure 9. - Concluded . 
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Figure /0.- History of events . Rigid model. 
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Figure //. -l/ariation of erperimentot mijcr wiitl di ' stance 1° max/mum 
occe/erofion increment for three conditions of test moc/e 1 and 
for four other gust -tunnel mode/s. 
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Figure 12 .- Comparison of calculated and experimental values of 


maximm ^ing- tip de flee t ion incre ment- for chap/p /p gradient distance, 
maximum deceleration incitement 
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(a) Sharp-edge gust with 2.4--ch ord gradient distance . Sample flight 2 . 
Figure 13. — Comparison. of experimental and theoretical reactions to a gust. F-I3.5cps. 
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(c) Gust with 21. 2-chord gradient distance. Sample flight 2. 
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Horizontal distance from front edge of gust tunnel, chords 
(a) 5harp -edge gust with 1.8 -chord gradient distance.Sa/nple flights. 
Figure 14- . — Comparison of experimental and theoretical reactions to a. gust, f = 2.6.1 cps. 
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( c ) Gust u/ith 22.4-chord gradient distance. Sample flight I. 
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Figure 15 .— Comparison of experimental and theoretical reactionstoa gust. Rigid wings* 
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(b) Rlode I A , condition £ . 

F/gure /6 - Continued . 
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(d ) Mode/ S , condition 2 

f/gure /6 . -Conf/nued. 
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Fig. 16e 
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(e) Hoc/e/ C , cond/t/on / . 

F/gure /6 . - Continued . 
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(f) Mode/ D } condition / . 

F/gure /6 -CcrtF/yue cf . 
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Fig. 16g 



(g) Mode/ D } condition 2 . 

F/gure /6 ~ Conc/uded . 
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F/gure 20 - Effect of change of w/ng frequency 
due to add /t /on of mass, node/ D , cond/t/ons 
/ and 2. 
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Figure 2L - Calculated history of reactions in a single gust. Mode/ C, condition I . 
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Fig. 22 





(o)H, =4.26 chords; b =4.27 . (b) H,= 10.25 chords-, b=i.26 . (c) H, =19.53 chords; b *0.459. 


Figure 22 Calculated history of reactions in a sing/e gust. Mode/ D, condition / . 
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Figure 23 . - Calculated history of reactions m a single gust. Model D, condition 2 . 
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Figure 24. ~ Calculated 


history of reactions in a repealed past. Mode/ C } condition / $ H t =fi 2 =9.99 chords 


tt 3 - 37 \ 46 chords ; U 2 - -U, 
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Figure 25.- Calculated history of reactions m a repeated gust. Mode! D, conditio/? /; H,~H 2 = /0.25 chords , 

H 3 =32.48 chords ; Uz m — U, . national advisory 
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Figure 26.-Colculated history 


of reactions 


in a repeated gust. Model D, condition 2$ H t - H 2 -10.08 chords, 

H, - 33.83 chords ; U z - -U, . 
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